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Abstract 



CN I A pedagogical presentation of integrable models with special reference to the 

Toda lattice hierarchy has been attempted. The example of the KdV equation has 
^ . been studied in detail, beginning with the infinite conserved quantities and going on 

On I to the Lax formalism for the same. We then go on to symplectic manifolds for which 

we construct the Lax operator. This formalism is applied to Toda Lattice systems. 

V^ . The Zakharov Shabat formalism aimed at encompassing all integrable models is 

^P ', also covered after which the zero curvature condition and its fallout are discussed. 

We then take up Toda Field Theories and their connection to W algebras via the 
Hamiltonian reduction of the WZNW model. Finally, we dwell on the connection 
between four dimensional Yang Mills theories and the KdV equation along with a 
Q^' generalization to supersymmetry. 
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1 Introduction: Non-Linear Equations 

Linear partial differential equations, in particular the Schroedinger, Klein-Gordon and 
Dirac equations, have been known in field theory over a long time, and have been used in 
many different problems with great success. Non-linear equations, i.e., equations where 
the potential term is non-linear in the field {S), have been known for some time as well. 
These equations and their solutions are the topic of the present Article. 

The earliest non-linear wave equations known in physics were the Liouville and Sine- 
Gordon equations. The Liouville equation arose in the context of a search for a manifold 
with constant curvature. Pictorially, such parametrizations may be likened to covering a 
surface with a fishing net. Since the knots on the fishing net do not move, the arc length is 
constant. The threads in the net correspond to a local coordinate system on the surface. 

The Liouville manifolds may be reparametrized locally so as to have a metric of the 
form: 



exp p , 

so as to be conformally equivalent to a fiat space metric. The study of such manifolds 
with constant curvature led J. Liouville [1] to the equation known by his name: 

expp (1.2) 



dxdy 

X and y being local orthogonal coordinates. Interest in this equation was renewed in the 
70's and 80's due to its appearance in string theories [2,3,4]. 

The Sine-Gordon equation, named after a pun on the Klein-Gordon equation, is an 
equation for the angle u between two coordinate lines when the total curvature is constant 
and negative. This equation first appeared in the work of Enneper in 1870, and has the 
form: 

— — = sin(a;) (1.3) 

oxoy 

where x and y are coordinates in a system with constant arc length. 

The Sine-Gordon equation has some interesting solutions known as solitons and 
breathers. A soliton satisfies three conditions. First, a single soliton must have con- 
stant shape and velocity. Secondly, it must be localized, and its derivative must vanish 
at infinity. Thirdly, if two solutions collide, they should survive the collision with their 
shapes unchanged. 

Principally, there are two types of solitons, one which increases by a fixed amount 
(say 2it), and is called a 'kink'; the other which decreases by the same amount, and is 
called an 'anti-kink'. 

A breather is a localized solution that varies periodically, and could be considered as 
a permanently bound system of a kink and anti-kink. 

An interesting property of the Sine-Gordon equation is that its solutions can be 
mapped into others through the Baecklund transformation [5], and can thus be used 
to create new solutions from known solutions. It is however impossible to generate a 
complete set of solutions from one original solution, via the Baecklund transformation [5] . 

A third non-linear equation which we shall study in some detail, was discovered in 
1895 by D.J.Korteweg and G. de Vries [6], while trying to describe the motion of water- 



waves in a canal. It has the form: 

ut - 6uu^ + u^^^ = (1.4) 

and is also known as the KdV equation. It has been extensively studied, and many of 
the properties of non-linear wave equations that are known today, were discovered in 
connection with its solution. This equation was solved by Gardner, Greene, Kruskal, and 
Miura in 1967 [7-13]. Along with N.J. Zabusky and C. H. Su, they also found many 
interesting properties of the same. One of these is that the KdV equation has an infinite 
number of conservation laws, and that the conserved quantities of each of these laws can 
be used as a Hamiltonian for an integrable system. This collection of Hamiltonians is 
called the KdV hierarchy. 

There exists a theorem of classical mechanics, which states that if a Hamiltonian 
system with 2n degrees of freedom has n functionally independent conserved quantities 
such that the Poisson bracket of any two of them vanishes, i.e., the integrals of motion are 
in 'involution', the system is completely integrable. It is clear that solutions of systems 
with an infinite number of conserved quantities must be infinitely restricted. A soliton 
is precisely such a solution: it is a localized wave which retains its shape even after 
collisions. Intuitively, it is clear that for this to happen, there must be an infinite number 
of conservation laws, and therefore an infinite number of conserved quantities. The terms 
'integrable models' and 'solitons' are often used synonymously. 

A system of coupled equations of motion describing a 1-dimensional crystal with non- 
linear coupling between nearest neighbour atoms, was introduced by M.Toda [14] in 1967. 
The equations of motion are 

m-^ = a[2e-''" - e^''-' - e~''"+'] (1.5) 

where r„ = Un+i — Un, and ■u„(t) is the longitudinal displacement of the n-th atom with 
mass m from its equilibrium position, a being a constant. These models admit soliton 
solutions which have been studied experimentally on an electrical network by K. Hitota 
and K. Suzuki [15]. In the continuunm hmit, these equations reduce to the KdV equation 
[5]. 

We see that models with exponential interactions are a source of non-linear equations, 
the Liouville and Sine-Gordon equations being examples. The Liouville equation could 
be generalized to include a mass term: 

^'"^ +m'<P = e^ (1.6) 



dxdy 

while the Sine-Gordon equation could be generalized to the " Sinh-Gordon" equation with 
the replacement uo ^ iu. Thus 

+ m ix) = smhu; (l-7j 



dxdy 
We also have the Toda Field Theory equations 



oxoy 



Here kij is the Cartan matrix for some complex Lie Algebra. The simplest of these field 
theories is the A.^. Toda field theory, and it includes the Liouville field theory for the 
special case r = 1. There exist generalizations of the Toda equations called " Afiine Toda 
Equations", and have an extra term on the RHS, taking the form: 



dxdy 



^e'^^.'Pj ^^R.e'^o'Pj (1.9) 



Here K is an affine Cartan matrix, and Ri the right null vector for this matrix when Rq 
is normalized to unity. 

These models include the Sinh-Gordon equation as a special case. Both the Toda 
and Afiine Toda field theories have an infinite number of conserved quantities [16]. They 
admit soliton solutions with an imaginary 0, [17]. Both models have been formally solved 
by Leznov and Saveliev [18]. 

The Toda field theories can be obtained from the Toda Lattice by setting 

i^i = {(pi - 0j-i) - (0i+i - (pi) (1-10) 

whence 

for SU{n + 1), showing that the space-independent solutions of (1.11) satisfy (1.5). 

Since the Toda field theories are the 7 = limits of the Afiine Toda field theories, 
they could be used to classify 2-dimensional models with a second order phase transition, 
with the Toda field theory describing the model at the critical point where it has to be 
conformally invariant [19]. Hence the great interest in (Afiine) Toda field theories. How- 
ever the precise connection is still unclear. Central charges and critical exponents have 
been calculated and compared. One hopes that the Afiine Toda field theories are per- 
turbations that correspond to the physical model away from the critical point. However, 
more explicit connections are yet to be found. 

The method originally used for solving non-linear equations, and especially the KdV 
equation, was the inverse scattering method originated by Gelfand and Levitan [20] . This 
involved looking for a linear equation related to the original non-linear equation, and 
studying the evolution of the latter. In 1968, P.Lax provided this method within a solid 
theoretical framework [21]. The Lax equation is 

Lt + [L,M]=0 (1.12) 

where L and M are operators satisfying 

Li) = \ij- (1.13) 

and 

^t = Mip (1.14) 

where A is a scalar, and ip a solution of a linear equation which is just the Schroedinger 
equation for the KdV case ! The Lax equation was generalized to the form of a zero 
curvature condition which facilitates greatly the form of the transition matrix from the 
initial to the final state. 



In what follows, we attempt to give a pedagogical presentation of Integrable Systems 
with special emphasis on the KdV and Toda systems. After an introduction to the KdV 
equation and its properties, we show how an infinite number of conserved quantities arise 
via the Muira [8] transformation, while detailed calculations are referred to ref.[22]. We 
then dwell on solutions of the KdV equation via the inverse scattering method and the 
Lax formalism [21], after which we obtain the Lax operator for symplectic manifolds, 
using the Toda Lattice as an example. The group structure of the Toda equations for 
SU{N) is also studied. The Lax transformation was later generalized by Zakharov and 
Shabat [23] to a first order formalism which was used by Ablowitz, Kamp, Newell and 
Segur (AKNS) [24], for a unified description of other integrable models. The essential 
features of this approach are also discussed. A fall-out of the above is the 'zero curvature 
condition' that facilitates the transition to the quantum case. However, the treatment we 
follow is strictly classical. 

Next we take up the Toda field theories, and after reporting briefiy the connection 
with conformal invariance, dwell on the Hamiltonian reduction of the WZNW model to 
the Toda field theory, which in effect transforms an affine Lie Algebra to a W-Algebra. 
(Most calculational details are skipped, but may be found in the literature [25]). Finally 
we refer to the interesting connection between the 4D self-dual Yang-Mills theory and 2D 
Integrable models, and the generalization to SuperSymmetry. 

The material is presented as follows. In Sect. 2, we introduce the KdV equation and 
its conserved quantities. In Sect. 3, solutions of non-linear equations are taken up, in 
particular the inverse scattering method and the Lax formalism. In Sect. 4, we digress 
to Symplectic Manifolds and construct conserved quantities for these manifolds. Sect. 5 
applies the above framework to the Toda Lattice where the group structure of the Toda 
equations is also discussed. In Sect. 6, we take up the unifying first order formalism of 
Zakharov and Shabat [23], continuing in Sect. 7 to the zero curvature formalism and its 
ramifications. In Sect. 8, we take up Conformal Invariance, and introduce Toda Field 
Theories which are constructed independently of the Toda Lattice. In Sect. 9, we carry 
out the Hamiltonian reduction of the WZNW model to Toda Field Theories. Finally 
in Sect. 10, we take up the connection of Toda Field Theories with Self-dual Yang-Mills 
models. Sect. 11 contains some concluding remarks. 

2 The KdV Equation 

The KdV equation was formulated to explain the solitary water waves observed by J.Scott 
Russell in the Edinburgh Glassgow canal. It is a non-linear equation in one space and 
one time dimension and possesses soliton solutions. Of this, however, nothing was known 
at the time of its formation. 

The KdV equation after an initial scaling takes the form 

du du d^u , , 

dt dx dx^ 

This equation is Galilean invariant, but not Lorentz invariant. It can be derived from 
the Hamiltonian 



where the u{x) satisfy the Poisson bracket relations 

[u{x),u{y)] = dj{x - y) (2.3) 

However, the Lagrangian from which it can be derived, is non local: 

1 |-+oo .du(y) f , .u^ l.du.o, , ,. 

LKdv = ^J dxdyu[x)t[x - y)—- J ^^~a ~ 2^^' ^ *- '' 

where 

e{x-y) = 0{x-y)-]^, (2.5) 

9 being the step function. Ergo, one cannot write down a local Lagrangian whose Euler- 
Lagrange equations yield the KdV equation. 

Solutions of the KdV equation can be shown to be soliton solutions which travel 
without any change of shape. It is the non-linear term which is responsible for the above 
property. 

What is most interesting about the KdV equation is that it admits of an infinite 
number of conserved quantities as was shown by Miura [8] . This procedure is explained 
below. 

The KdV equation is related to another equation called the modified KdV {MKdV) 
equation, viz., 

dv 2^"^ 9^"^ 
dt dx dx^ 

where v is related to u in the KdV equation through the Riccati transformation 

u = v^ + iVQ^ (2.7) 

ox 

The MKdV equation is however not Galiliean invariant. Under the transformation 

3t 3 ve \^ 

it reduces to 

dtv = (^ + v)d.v + dlv = 9.[^ + ^ + dlv] (2.9) 

This yields a solution of the KdV equation through the transformation 

M = eV/6 + t; + ieS^t; (2.10) 

The second form of (2.9) is in the nature of a continuity equation, so that we can identify 

/+00 
dxv{x{t)) (2.11) 

-00 

as the conserved quantities, v can be inverted in terms of u as 

00 
V = Y.t''VnHx,t)) (2.12) 



and this yields Vn{u{x,t)) as the conserved densities, since each power of e must indepen- 
dently satisfy a continuity equation. That these are also in involution can also be checked, 
being explicitly shown by Das [22]. Some of the conserved quantities are 

u u 

vi = -id^ui] V2 = — - - dlu] V3. = id^— + dlu] (2.13) 
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3 The Lax Framework 

Linear Hamiltonian systems with fixed initial value problems can be solved using the 
Laplace or Fourier transformations. Such methods are inapplicable for the nonlinear 
equations and new methods must be found. Gardner, Green, Krushal and Miura [9] 
managed to solve the initial value problem for the KdV equation in a very ingenious way. 
In subsequent years, this method has become the standard method for solving non-linear 
systems and goes by the name of inverse scattering theory [20,21]. This method is outlined 
in Fig 1. 

The initial value for the partial differential equation is used as the potential in a 1- 
dimensional scattering problem for a linear equation, e.g. the Schroedinger equation. One 
then finds the so called scattering data, i.e. discrete spectrum, normalization constants, 
reflection constants (as a function of the wave number) for this scattering problem. Using 
the partial differential equation [pde] evaluated for \x\ asymptotically large, (and hence 
the "pde becomes a linear equation because the potential is assumed to vanish at spatial 
infinity), the values of the scattering data can be found for all later times. Finally, the 
scattering data allow one to reconstruct the potential, and hence the solution of the prfe 
for any later time. 

One would intuitively like a better understanding of the origin and relevance of the 
linear Schroedinger equation. One way to see this is through a generalized Riccati relation 
of the form: 

.2 , ,.la-- (3_l^ 



U 



6A 



ox 



so that the KdV relation (2.1) reduces to 



dv 

m 



6A) — 
dx 



dx^ 







(3.2) 



As mentioned earlier, a solution of the MKdV equation yields a solution of the KdV 
equation through the Riccati relation. The simplest way to attempt an inversion of the 



Riccati relation is to linearize it. To that end we define 

V = iyQipx/il} (3.3) 

so that (3.1) takes the form 

M + 6A = -6^,,/^, (3.4) 

or equivalently, 

u 

^.. + (- + A)V = (3.5) 



which is, in fact, the time-independent Schroedinger equation. There exists however a 
more formal theory due to Lax [21], which we now elaborate. 

Given a linear equation described by a time-independent Hamiltonian if, and an 
operator A whose expectation values are time independent, A{t) is unitarily equivalent 
to A(0): 

\j\t)A{t)U{t) = A{Q) (3.6) 

where U{t) is the time-evolution operator with the form 

U{t) = exp[-iHt] (3.7) 

Differentiating (3.6) gives 

8A 

U\t){—-t[A,H])Uit) = (3.8) 

which implies 

— = z[A,i7] (3.9) 

Thus for the expectation value of A{t) to be time independent, the standard time evolution 
relation (3.8) must be satisfied. Further, from eq.(3.7) follows the relation 

^^ = -tHU{t) = BU{t) (3.10) 

where 

B = ~iH (3.11) 

is an anti-Hermitian operator. 

This argument is mimicked in the case of a non-linear evolution equation. Let 

L{u{x,t)) = L{t) (3.12) 

denote the linear operator we seek. We assume it to be Hermitian, and to have eigen- values 
independent oft. For this to be true, one must have u^{t)L{t)u{t)=L{Q). Differentiating 
both sides w.r.t. t, we obtain 

^^ui),uHi)'-^im.imm^ ,3.13) 

Unlike the linear case, we do not know the form of U{t). However, U is unitary, so 

dt ^ ' dt ^ ' 



Thus we can write 

^ = Bmit) (3.15) 

where anti-hermiticity must be imposed on B. Substitution in (3.12), and a httle simph- 
fication, yields 

^ = [B{t),L{t)] (3.16) 

which is similar to (3.8), except for the fact that we do not yet know the form of B. 
However, let us assume that L{t) is linear in u{x,t). Consequently, the LHS of (3.14) 
is a multiplicative operator, proportional to the time evolution operator of u{x,t). This 
would ensure that the eigen- values A of L(t) would be time-independent, i.e., 

L{t)ilj{t) = -\i){t) (3.17) 

Further, ip{t) must be unitarily related to its value at t = 0, i.e., 

V^(t) = f/(t)V^(0), (3.18) 

and its evolution w.r.t. time would take the form 

dtl){t) dUit) 



dt dt 



-^PiO) = B{t)m (3.19) 



The operators L{t) and B{t), when they exist, are known as the Lax pair, corresponding 
to a given non-linear evolution equation, and play a fundamental role in determining the 
solution. For the KdV equation, L[t) is obtained from the linear form of the Schroedinger 
equation 

L{t)=D^ + ^u{x,ty, D = ^. (3.20) 

By trial and error, B{t) can be chosen so that (3.15) is satisfied, and a possible solution 
is 

B{t) =4:D^ + -{Du + uD) (3.21) 

The solution for ip w.r.t. t follows from (3.18) and (3.20) to be 

1 

i^t = '^i'xxx + ^U^i' + U^Ja: + COUSt.^p (3.22) 

which yields, using the Schroedinger equation (3.5): 

ipt + -Uxi^ - -ui^x + 4A^^ = const.ij (3.23) 

6 

A.Lenard [26], in an unpublished report, further displayed the relation between the 
Schroedinger equation and the KdV relation by elegantly deriving the latter from the for- 
mer, using only the assumption that the spectral parameter A in (3.4) is time-independent. 

The KdV equation exhibits also a fascinating symmetry, i.e., that of the group 
SL{2,R). Consider a group element 

g = exp[te-Ta] (3.24) 



where Ta is a generator of SL{2, R), and define 

A^ = g-'d,g (3.25) 

Then the KdV equation follows from the fact that the Maurer-Cartan equation 

d^A,-d,A^-[A^,A,]=0 (3.26) 

is satisfied for a special for a special choice of gauge, e.g., 

Al = -v^; (A < 0); Al = 6; Aj = -^uix, t); A^ = A(m(x, t)) (3.27) 

4 Lax Formalism On Symplectic Manifolds 

In this Section, we conclude the above study of the KdV equation with a with a short 
discussion on symplectic geometry, which is directly relevant for application to the Toda 
Lattice. 

A symplectic manifold is one with a preferred 2-form /^j, which is non-degenerate and 
closed. The phase space of an integrable model corresponds to a very special symplectic 
manifold, since it possesses a dual Poisson bracket structure. We assume that there exist 
two distinct 2-forms which are both non-degenerate and closed. One way of expressing 
the existence of two distinct symplectic structures is to require that the same dynamical 
equation be described by two distinct first order Lagrangians Lq and L, where 

Lo = 9^^\y)r - Hoiy); (4.1) 

L = 9,iy)r-Hiy) (4.2) 

where 

r = ^; [/x=l,2,...2iV] (4.3) 

The Euler-Lagrangian equations following from (4.1-2) are 

Uuiy)y'' = d.Hoiy) (4.4) 

F.My" = d,H{y) (4.5) 

where 

f,. = d,e^\y)-d,ef{y) (4.6) 

F,. = dMy)-d.O,{y) (4.7) 

It is easy to see that the two forms / and F are closed, where 

f = \f,.dy^^dy^ (4.8) 

F = \F^.dy^ A dy"" (4.9) 

since f^y and F^^ satisfy the Bianchi identities 

5a/m. + -^^/.A + 5,/a^ = 0; (4.10) 

10 



and 

dxF^, + d^F^x + d^Fx^ = 0. (4.11) 

Besides, they must also be non-degenerate since (4.4) and (4.5) describe the same dynam- 
ical system. Let their universes be f'^'^ and F'^", i.e., 

f,ur^ = F,,F^^ = 5l (4.12) 

so that (4.4-5) take the forms 

y^ = rd,H,{y) (4.13) 

r = F^^d,H{y). (4.14) 

We can also construct a nontrivial (1,1) tensor S*!^ as 

s; = F,Mf"'iy)- (4.15) 

Consistency of (4.4) and (4.5) further requires that 

d,d,Ho{y) - d,d,Ho{y) = (4.16) 

so that after a little algebra, one can show that 



dt 
where 



-U'h^ + U^fx, (4.17) 



U; = d.y'' = d,[r'dxHo{y)] = d,[F'''d,H{y)] (4.18) 

with a corresponding relation for F'^'^, i.e., 

^^^ = -U^F,^ + U^F,, (4.19) 

involving the same [/-tensor. The corresponding equations for the inverses f^'^ and F^'^ 
follow from (4.17) and (4.19), and have the forms 

JfflU 



We can finally show that 



which in matrix notation 



dt 


-J ^A - 


-J ^A 


dF^""" 
dt 


= F^'^U^ - 


- F'^^f/j; 


ds;: 

dt 


= s^ui - 


- u^si 




tHS,., 



(4.21) 

(4.22) 
(4.23) 



can be recognized as a Lax equation (3.15), thus providing a Lax representation of the 
dynamical equations (4.13) and (4.14). One important consequence of (4.23) is that the 
set of quantities 

Kn = -TrS"" (4.24) 



n 
11 



fix, y) = Tr3 + oid.u + ud,)6{x - y) (4.30) 



and 

Ko = ln\ detS \ (4.25) 

can be shown to be invariants since 

^ = ^^[^(^)§] = Tr[P{S)[S,U]] = (4.26) 

P{S) is a polynomial in S. That these are in involution can easily be checked, as done 
explicitly in ref. [22]. Applied to the KdV equation, the two Poisson structures of that 
equation are given by the correspondence: 

F^^ -^ D; (4.27) 

Going to the coordinate bases we have 

F{x,y)=<y\D\x>=dxS{x-y) (4.29) 

d^ 1 

so that 

F-\x, y) = e{x - y) = e{x - y) -]^ (4.31) 

However f~^{x — y) cannot be expressed in a closed form. The Lax operator S takes the 
form 

S = D^ + -u+\{Du)D-\ (4.32) 

and with a little algebra, (4.23) can be shown to be reduced to the KdV equation, with 
consequently an infinite # of conserved quantities. This is described in detail in ref. [22]. 

5 The Toda Lattice 

The model of the KdV equation that has been studied so far is a continuum model. A 
finite dimensional system with a finite # of degrees of freedom is simpler to study. The 
Toda Lattice is such a system to which the symplectic approach of the above Section 
is especially applicable. We now study the Toda Lattice and its integrability from a 
symplectic point of view, following it up with a group theoretical treatment. 

The Toda Lattice describes the motion of N point masses on the line, under the infiu- 
ence of an exponential interaction. The Hamiltonian equations in terms of the canonical 
coordinates Qi and momenta P^ are given by 

(5.1) 



g, = P,; (2 = 1,2,. ...AT); 






p. = Q-(Q,-Qj-i) _Q-'(Qi+i-Qj). 


{j = 2,3,.N- 


-1); 


Pi = -e-(«2-«^); P^ = e-(«-- 


-Qn-i)_ 





The equations can be cast into a more symmetrical form by enlarging the system to 
(A^ + 2) point masses, with end points at spatial infinity. In that case, the Hamiltonian 
equations take the form : 

4 = p.; (i = l,2,....Ar);Pi = e-('3"-'3")-e-('3"-Qi. (5.2) 

12 



We can choose 

y' = Q^■, y''-'' = P.; (^ = i,2,..iv). (5.3) 

Applying the geometrical method of the previous Section, two choices of the Lagrangian 
are as follows: 

Lo = J:kiP^Q^ - Q^Pi) " T^Pf + e-(«»+-«»)]; (5.4) 

TV -, 

L = T\o{P- + e-(«'+-«>))Q. + 7r,(P)P,] - H{Q, P) (5.5) 

1=1 

where 

1 ^ 

TV p3 

i/(g, P) = ^[^ + (P, + p^^,)e-(Q-+i-Q^)] (5.7) 

/^i, turns out to have the canonical Poisson bracket structure 

so that 

r = i^j o) (^-^^ 



F^y can be shown to have the form [22] 



B e 
where 



F,.= (^ ^) (5.10) 



Ai^. = 5^_^^_^.e-(Q'+i-Q')-5,,j.^ie-(«^+i-Q^) (5.11) 



The (1,1) tensor S^ thus takes the form 



s; = (^, p) (5-12) 



and the conserved quantities are 



TV 

TrS = 2TrB = 2 ^ P^; (5. 13) 

1=1 

1 ^ P2 

-TrS^ = Tr[2p2 - {Ae + eA)] = Y,[— + e-^^'+i-'^.)] = ^^(q^ p). (5,14) 

2 j=i 2 

1 Af p3 

TrS^ = J2[^ + {P, + p^^^)e-W«+-«-)] ^ iJ(g, P) (5.15) 

13 



The Lax representation (4.23) for the Toda equation takes the form of the following matrix 

equations 

d A 

^ = -lB.Dy, (5^16) 

^^^A-De^lle.D] (5.17) 

which reduce to the Toda equations Qi = Pi and Pi = e~('3«-Qi-i)_g-(Qi+i-Qi) respectively. 

5.1 Group Structure of Toda Equations 

Eq.(5.1) can be differentiated and put in the form 

(5.18) 

-(Qi + l—Qi) 



Qi - 

Qn -- 


= _g-(Q2-Qi) 

= p_e-(Q.-0.-i)_g 

= Pn = e-(«^-«^-i) 




N N 
1=1 i=l 



It is easily checked that 

N N 

" ' (5.19) 

i.e., the total momentum is conserved, and therefore the centre of mass motion can be 
separated and the dynamics of the system expressed in terms of (A^ — 1) coordinates and 
momenta. Defining 

qa = Qa+i-Qa; a = 1, 2...iV - 1, (5.20) 

the second order equations satisfied by the g^'s can be written as 

gi = 2e-'^' - e-'^' (5.21) 

Qa = -e-"^-' + 26-"^ -€'"''+'■, a = l,..N-l 
(JN = -e-*?^-! + 2e-«'^ 

which can be compactly written as 

N-l 

qa = E KabC-'' (5.22) 

b=l 

Kab being the Cartan matrix for SU{N). Eq.(5.22) generalizes for the other Lie Algebras 
as well. 

The Lagrangian giving rise to the above Euler-Lagrangian equations can be written 
as 

N N ^ N 

L = J2J2 7;9aKab-Hb - E e-^'' (5.23) 

a=l b=l a=l 

Kab~^ being the inverse of the Cartan matrix. The momenta conjugate to qa are defined 
as 

QT N-l 

Pa = ^=J2 Kab~'qb (5.24) 

and it is easily checked that 

{qa,Pb} = 5ab (5.25) 
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so that {qa,Pa} constitute a canonical coordinate system. 

That the group structure entering above is not just accidental, can be seen by defining 
the following Lax operators: 

1 ^ 
S=- Y.[PaHa + {Ea + E^a)e-'^^/\ (5.26) 

^ a=l 

1 N-l 

t/ = --Ee"'"^'[^a-^-] (5-27) 

^ a=l 



where Ha and Ea are the generators of SU{N) in the Chevally basis. 

The Lax equation (4.23) can be seen to be satisfied, since ^ — [S, U] reduces to 

1 N-l N-l 

- Y. HaKab~'[qb-T.^bce-'^] 

^ a,6=l c=l 

which is zero by virtue of the Toda equations (5.22). Hence the quantities 

K^ = -TrS" (5.28) 

n 

must be conserved under the flow of the Toda equations. Since S belongs to the SU{N) 
algebra, the number of independent conserved quantities can equal {N — 1), which is the 
rank of SU{N). The total number of conserved quantities is thus A^, if we add the total 
momentum. It can be shown that these are also in involution [22]. This treatment is due 
to Leznov and Saveliev [18]. 

6 Zakharov-Shabat Formalism 

So far we have only studied two integrable models, viz., the continuum KdV and the flnite 
dimensional Toda Lattice. In trying to understand the non-linear Schroedinger equation 
which is also integrable, Zakharov and Shabat [23] obtained a description which was later 
generalized by AKNS [24] to describe various other integrable models. This approach uses 
a Lax operator which is flrst order in the derivative d^, in contrast to the second order 
formalism in eq.(3.19). Besides describing various integrable models in a unifled manner, 
this approach has the additional advantage that the inverse scattering method generalizes 
readily to the quantum case. In what follows, we describe the flrst order formulation of 
the Lax operator, and elucidate the essential features of this approach. It is easily checked 
that if 

L(t)^(t) = -A^(t); (6.1) 

dtm = [B{t),L{t)], (6.2) 

where 

dij{t) 

then 



Bit)m, (6.3) 



^-^ 
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We can invert the argument to identify the Lax pair in the following way. Namely, if 
mm = -A^(t); ^ = B{t)i,{t),with^ = 

, i.e., if the compatibility condition of (6.1) and (6.3) yield the system under study, then 
L(t) and B{t) can be identified as the Lax pair of the system. We would like L{t) to be 
linear in d^- Using the analogy between the Klein-Gordon and Dirac equations, we define 
a two-component column matrix 

and generalize the two equations (6.1) and (6.3) to first order matrix equations as 

— = (go-+ + ro-_ipo-3)0; (6.6) 

^ = {Pa+Qa+Ras)(i) 

where (Tj- and a^ are the Pauli spin matrices. The dynamical variables q{x,t) and r{x,t) 
do not depend on the spectral parameter p which is assumed to be independent of x 
and t. The coefficient functions P, Q and R on the other hand, do depend on p, and 
are functionals of q and r. Demanding that the partial derivatives of w.r.t. x and t 
commute, we obtain the compatibility conditions to be 

-g = <lQ- rP; (6.7) 

|.^_2.«_2,,0; ,6.8) 

| = f + 2,fi + 2.pP (6.9) 

i.e., if (6.7-9) describe the non-linear evolution of a system, then (6.6) describes the Lax 
pair appropriate for such a system. Explicitly 

L = dx — qa^ — ra^ (6.10) 

B = Pa+ + Qa^ + Ras (6.11) 

so that (6.2) is satisfied. 

The choice of r = 6 yields the KdV equation, and the choice r — q = — if /^/6, the 
MKdV equation. The choice q = vk^p* and r = vk^, k being an arbitrary constant 
parameter, yields the non-linear Schroedinger equation: 

tdt^ = -^,, + 2k\^fij (6.12) 



and the choice r = — g = ^Ux with 



P = Q = -— sincj 
Ap 



yields the sine-Gordon equation. 
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The operator (L + A) in (6.1) can be rewritten as f (x, t, A) + dx, where 

V = ~q(J+ — rcr^ + ipo^ (6.13) 

If one knows the solution of the associated Schroedinger equation at some other point 
{x,t) by multiplying the solution by a hermitian matrix T{x,y,t, A), i.e., 

^ljix,t,\) = Tix,y,t,\)^ljiy,t,\) (6.14) 

where T(x, y, t, A) is a solution of 

dxT{x,y,t,\) = -{q(J+ -ra_ + ipa3)T{x,y,t, \) (6.15) 

with the initial condition T{x, x, t, A) = /. 

7 The Zero Curvature Condition 

The Lax condition (6.2) can be written as 

[dt-B,L]=0 (7.1) 

Using 

L = dx- A{x) (7.2) 

we obtain the form 

m-B),idx-A)]=0 (7.3) 

which is like a zero-curvature condition for 

Foi = [(9o-Ao),(9i-A)] (7.4) 

with the identification 

Ao = -B{x, py,A, = -A{x, p) (7.5) 

The importance of the zero curvature condition stems from the fact that (6.6) may be 
solved, using 

iIj{x) =T{x,y,p)ilj{y) (7.6) 

where the transformation 

T{x,y,p) = P,exp[- f A^{z)dz\ (7.7) 

where P^ denotes path ordering. 

It is easy to see that T{x,y,p) translates solutions of the problem along the x-axis 
for a fixed time, i.e., 

T{x,y,p)T{y,z,p) = T{x,z,p); (7.8) 

T'^{x,y,p) = T{y,x,p); (7.9) 

T{x,x,p) = l (7.10) 

Setting 

Ur{x2,t2;xuti) = P^exp[- H'' A^dx^"] (7.11) 

J Xl,tl 
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and taking the product of two such exponents, it is easy to see that 

Uri{x2,t2; xi,ti)Ur2{xi,ti; X2,t2) = exp[-- j da^'^F^y], (7.12) 

using the Baker- Campbell-Hausdorff formula and the Stokes theorem, the integration 
being done over the area enclosed by the closed path ri + r2. As the curvature F^j_i, 
vanishes, 

Uri{x2,t2; Xi,ti)Ur2{Xi,ti; X2,t2) = 1 (7.13) 

and so 

U^^ix2,t2;xi,ti) = UriXi,ti;X2,t2) (7.14) 

so that 

Ur^{x2,t2;Xi,ti) = Ur^{x2,t2]Xi,ti) (7.15) 

ergo, U is independent of the path taken. For a closed path, f/(a;, t; x, t) = 1. Hence path 
ordering drops out of the transition matrix T{x,y,p). 

Returning to the time evolution of the transition matrix, it can be shown that 

dtT{x, y, p) = [B{x, p),T{x, y, p)] (7.16) 

which is the form of a Lax equation, so that all quantities of the form 

ir„ = ^Tr[T{p)r; K, = \n\detT{p)] (7.17) 

are conserved. We thus have an infinite number of conserved quantities when the zero 
curvature conditions are fulfilled. That this holds also for Toda Field Theories was shown 
by Olive and Turok [16]. 

8 From Conformal Invariance To Toda Field Theory 

That the KdV equation has a hidden conformal symmetry can be seen by making a 
Fourier expansion with Fourier coefficients 

1 z'^''" dx 

nn = —J tx(x)e-^"V^"0 (8.1) 

4 Jo 2,71 

It can be shown that the Poisson brackets of the «„ satisfy the Virasoro Algebra (up to 
trivial factors), i.e., 

- 2i7r{M„, Um} = -(n - m)u„,+n + ^^(^^ ~ l)^n+m (8.2) 

Higher order terms in the KdV hierarchy have a hidden u symmetry. 

We now digress to take a look at Toda Field Theories. These are essentially the 
only class of integrable, interacting, conformally invariant field theories in two space-time 
dimensions. To see this, we start with the generic action 

S = J[^^^<|)^^y^ - V{<j),)]d^Z (8.3) 



The trace of the naive conserved energy-momentum tensor becomes 

T^ = 2V. (8.4) 

As the trace of the energy-momentum tensor is required to vanish in a conformally 
invariant theory, it seems that if V 7^ 0, the theory is not conformaL However there is 
an ambiguity in the definition of the energy-momentum tensor. If we attempt to improve 
the naive energy-momentum tensor without violating the conservation property, we could 
choose 

V = ^M- + [^^.^u - r?^.5V(0.) (8.5) 

whence the trace of the modified energy-momentum tensor is 

e^^ = 2V + d+d-f (8.6) 

± being the light cone directions. If the second term is to cancel the first, we somehow 

need to get rid of the derivatives. This can be done, using the equations of motion. 

Without knowing the explicit equations of motion, the most general expression for /(0j) is 

Y^ Ci(j)i. Using the equations of motion resulting from varying the action , the tracelessness 

condition becomes 

dV 
2^ + Ec.7r-=0 (8.7) 

i d(f)i 

Eq.(8,7) is easily solved, with the result that the trace of the energy-momentum tensor 
vanishes if the potential is of the form 

¥{(/),) = Y.d,exp['£b,,<l),,], (8.8) 

j 

satisfying the requirement 

EcA, = -2 (8.9) 

i 

We choose bij to be related to the Cartan matrix of a simple Lie Algebra. The resulting 
field theories are called Toda Field Theories, and are described by the action 

1 2 

Sroda = j[-{d^<P,dy) - ^^exp(/? < a^\(j,>)](fx (8.10) 

where <, > is the scalar product in the root space, and takes its values in the root space 
of the simple Lie Algebra on hand. 

The equations of motion obtained from (8.10) are 

i3d>'d^(t)i + m^exp (^ Kij(t)j) = (8.11) 

Specializing for the SU{n) group, and setting m = f3 = 1, this becomes 

d+d.(P, = -exp {K,,(Pj) (8.12) 

With 00 = and 0i+i = 0, this reduces to 

d+d-(j)i = -exp (20i - 0i_i - 0i+i) (8.13) 
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Setting 

^j = {4>i - 4>i-i) - (0i+i - 4>i), 

after Mikhailov [27], we get the equation 

dt% - d.^ij, = -[2e^> - e^-i - e^>+^] (8.14) 

which is easily seen to be related to the Toda equations (5.20). One expects that the 
Toda Field Theories are integrable, and it turns out that they are indeed so (see ref.[8]). 
The calculation rests upon the existence of a zero curvature condition for certain group 
theoretical combinations of 0, which can be chosen as gauge fields. 

As mentioned earlier, the Toda Field Theories have been completely solved for simple 
g by Leonov and Saveliev [18]. They have also been solved for affine g by Olive and Turok 
[16]. 

Quantization of the Toda Field Theories is more problematic since the potential has 
no local minimum, the latter being attained at infinity, using the gauge group Ai. A lucid 
discussion of the problems encountered in the theory is given in ref. [29]. 

The central charge of the Toda theories can be constructed using free field technology, 
and is found to be [30] 

Zn 4:71 p 

r being the rank of the algebra, p being half the sum of the positive roots, and p^ its 
dual. Eq.(8.15) gives an indication that a quantum Toda theory with a strong coupling 
constant is equivalent to another Toda theory with a weak coupling constant, obtained 
by replacing j3 by Air/hf], and interchanging roots and "coroots". 

Incidentally, strong/weak coupling duality has recently become a subject of immense 
study in relation to string theories. 

It is possible to obtain the minimal models from the Toda Field Theories. For a 
particular value of (3, the central charges can be made to agree. However this is not 
enough. A complication arises from the fact that not all primary fields in the minimal 
models are actually present in the Toda theory. However, because of the duality in the 
theory, we can add another part of the potential with the coupling constant replaced by 
its dual; see Mansfield [31]. This modification is sufficient to give complete agreement. 

9 W- Algebras: Hamiltonian Reduction of WZNW 

Another fact which makes the conformally invariant Toda theories interesting is that to 
each such Toda theory, there corresponds a W^-algebra. The W^-algebras are an extension 
of the Virasoro algebra by adding primary fields primary fields of spin higher than Z, 
and were introduced by Zamolodchikov [32] as a pointer to conformal field theories with 
a larger overall symmetry. Zamolodchikov [32] investigated the case in which a primary 
field w{r) of weight 3 is added to the Virasoro algebra. In order for the algebra to be 
close, it had to be made 'non-linear', and hence lost its linear Lie Algebra character. 

Balog et al [33-35] showed that the Liouville and Toda Field Theories can be obtained 
as conformally reduced WZNW theories. This reduction can be viewed as a gauge 
procedure, and the Toda field theory can be obtained as the gauge invariant content of 
a gauged WZNW theory. The Liouville theory is obtained for the special case of the 
SL{2, R) gauge group. 
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The most powerful method of constructing W^-algebras is through the so-called quan- 
tum Drinfield-Sokolov reduction. In this, one starts with an afiine Lie Algebra, and 
reduces it by imposing some constraint on its generators. At the classical level, this 
procedure which leads to the so-called Gelfand-Dickey algebras [36], was pioneered by 
Drinfield and Sokolov [37]. 

It is thus clear that under the reduction that takes a WZNW field theory to a 
Toda field theory, the affine Lie Algebra that characterizes the WZNW theory reduces 
to a W^-algebra that is associated to a Toda field theory. This approach is also readily 
generalizable to the supersymmetric case where various new I^-superlagebras have been 
found as symmetry algebras of supersymmetric Toda field theories. We refer the interested 
reader to ref.[38] for further progress in this area. 

In what follows, we review the essential steps of the Lagrangian reduction of the 
WZNW model. The WZNW action for a non-compact group G in 2D Minkowski space- 
time is 

S{9) = -^ f d'p^^''Tr{g'%g){g~'d.g) + -^ f Tr{g-'dg)' (9.1) 

where B is the volume occupied by S'^. The left and right Afiine Kac- Moody [AKM] 
symmetries of this theory are generated by the Noether currents 

J(A) = KTr[\{d+g)g-']- J(A) = -KTr[\g-\d^g)] (9.2) 

where n = ^, and A is an element of the Lie Algebra g. The WZNW equations of 
motion are known to be equivalent to the current conservation 

a_J = (9+J = 0. (9.3) 

We now choose the following Gauss decomposition of an arbitrary element g=ABC, e.g., 

(9.4) 



A = 


= exp y^ x"Eo^ ; 

aeA+ 


B = 




C = 


= exp y^ y'^Ea ; 

aeA- 



where Cartan-Weyl root vectors Ea, Cartan subalgebra generators Ha = [Ea,E_a], and 
a set of positive (negative) roots A^ have been introduced with the following properties 

Kap = Oi{Hp) = - — ^; a,/3GA; | a^on^ |^ = 2; (9.5) 

I '^ I 

Tr{HaHfs) = -^Kafs = C^fs; (9.6) 

I <^ I 

Tr{EaEp) = -^5a,-(3; Tr[Ea, Hp] = 0. (9.7) 

I <^ I 

We also introduce the Polyakov-Wiegmann identity 

S{ABC) = S{A) + S{B) + S{C) + K I (fpTr[{A-^d-A)d+B)B-^ (9.8) 

+{B-^d^B){d+C)C-^ + {A-^d^A){B{d+C)C-^B-^)] 
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We now see, using eqs. (9.4-9.8), that the generahzed constraints 

J{Ea) = Kc"^; J{E^a = -Kc^; a G A+ (9.9) 

with some real numbers c^g whose values do not vanish only for primitive roots a G A, 
are enough to reduce the G- based WZNW theory to the Toda Field Theory defined by 
the Lagrangian 

Lroda = -^[\c^pd^rd^<pf' - E {u^e^'^^-^^')] (9.10) 



where 



(m2)« = I a l^c^c^ 



. Due to 0^2 7^ for the primitive roots, the constraint (9.7) can be re- written in terms 
of the Gauss decomposition (9.5-7) as follows: 

A-'d_A = B[J2hafc^E^]B-' (9.11) 



1 1 

E -I a fc^Eaexp[-Ko,p(l)% 

oGA 



id+C)C-' = 5-i[E J|a|'c?E„,]5 (9.12) 

aGA ^ 

= $:^|a|'c?i?_„exp[l/r,^0^]; 

oGA 

In the WZNW equations of motion, A and C occur only in the combinations given in 
(9.11-12), so that they can be eliminated in favour of B or 0". The remaining equation 
is just the Toda equation [25,34,35]: 

d+d.r + l\a fiu'^fexpi^K^pcP^] = 0; (9.13) 

(see also ref.[25] for details). 

As mentioned earlier, the Toda Field Theory possesses an extended symmetry rep- 
resented by a classical W^-algebra. These W^-algebras can be obtained as the quantum 
versions of the so-called Gelfand-Dickey algebras [36] known in the theory of KdV equa- 
tions. For instance, the Poisson bracket associated with the KdV equation in (8.2), results 
in the classical version of the Virasoro algebra which is the simplest W^-algebra. Moreover, 
the Lax representation of the KdV equation (3.15), defines the third order differential op- 
erator B = w^^\ The Fourier components of B, along with those of the KdV field, form 
the Gelfand-Dickey [36] algebra that generalizes to w^^^ in the quantum case. 

Now regarding the Toda theory as a constrained WZNW theory, the Hamiltonian 
structure can be obtained by a classical Drinfield-Sokolov reduction from the constrained 
phase space of the AKM algebra. In the Hamiltonian formalism, the AKM symmetry 
of the WZNW theory is represented by first class constraints. The W^-algebra of the 
Toda theory arises as the Poisson bracket algebra of gauge-invariant polynomials of the 
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constrained AKM currents and their derivatives. In what follows, we summarize the 
arguments supporting these statements. 

Let g{z, z) be the G- valued WZNW fields and J{z) the corresponding AKM currents 
having the form 

g{z, z) = g{z)g{z)- dg{z) = J{z)g{z) (9.14) 

Let dimg be the dimension of G; / its rank; k the level of the associated AKM algebra 
g\ g the dual Coxeter number of G\ p the half sum of the positive roots; and (3 the dual 
of p. 

The constrained WZNW theory is specified by (9.9). After a suitable choice of 
constants q, the currents J{z) can be decomposed as 

J{z) = K+jiz); K=Y,E_^- (9.15) 

i=l </>eA+ 

where {Ea^} are / simple roots of g. The maximal subgroup of G leaving this form of 
currents invariant, is the maximal nil-potent subgroup generated by E(f), (0 G A+), and 
implemented by the {dimg — l)/2 constrained AKM currents J'^{z). This allows us to 
interpret the constrained WZNW theory as the gauge theory in which all but / of the 
{dimg + /)/2 components of J are gauge components [33-35]. 

The current j{z) and the gauge transformations corresponding to E^ act on each 
column of the WZNW field g{z) separately, while each column contains only one gauge- 
invariant component e (of the highest weight), satisfying E^e = 0. The gauge degrees 
of freedom corresponding to the other elements of each column can be eliminated by a 
gauge fixing in favour of e. Because of (9.15), this leads to a linear pseudo-differential 
equation De = 0, where D is a polynomial pseudo-differential operator whose coefficients 
are gauge invariant polynomials in the currents J. This operator D can now be used to 
define a classical W^-algebra by choosing a Drinfield-Sobolov gauge in which one has 

where P's are the orders of / independent Casimir operators of g, and Fp generators with 
H weights (P — 1), so that the gauge-fixed current (9.16) has only one non- vanishing 
component in each of the / irreducible representations in a decomposition of the adjoint 
of g w.r.t. one of its sub-groups SL{2,R). The Poisson brackets between the different 
polynomials w^ define a classical W^-algebra. 

We close this Section by noting that Toda field theories also play an important role 
in the discussion of VT-gravity, where they arise as effective quantum theories [39,40] for 
the FT-gravity degrees of freedom in the conformal gauge. For a quantum version of the 
WZNW — > Toda conformal reduction, see [34, 41]. 

10 Self-Dual Y-M Theories: 2D Integrable Models 

The self-dual Yang-Mills (SDYM) theory appears to be a master theory for a whole variety 
of 2D integrable systems, as we are now going to explain. Though there is no general 
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proof, the statement can be checked on a case by case basis. The main point is that the 4D 
self-duahty condition admits of a zero curvature representation underlying a Hamiltonian 
description of SDYM descendents in lower dimensions. This makes it possible to apply 
the inverse scattering method for integration of the SDYM equations. Simultaneously, it 
explains the origin of gauge symmetries in integrable systems of the KdV type, since the 
SDYM theory in both gauge and conformally invariant in 4D. And last but not least, this 
connection provides us with a systematic way to associate the KdV type hierarchy with 
any simple Lie Algebra. 

SDYM solutions invariant by the action of a subgroup with two conformal generators 
satisfy a 2D differential equation, since eacd ID subgroup reduced the number of indepen- 
dent variables by one. This allows us to describe the invariant SDYM solutions in terms 
of a 2D integrable system. All known 2D integrable systems seem to be derivable this 
way, by appropriate truncations of a 4D self-dual gauge theory. This is true, in particular, 
for the KdV and non-linear Schroedinger equations, the Liouville and Toda equations, as 
well as other integrable in 2 and 3 dimensions. Our presentation in this Section is only 
illustrative; we give one explicit example of embedding of the KdV equation into the 4D 
SDYM theory [42], and a supersummetric generalization. 

Let x"' = (x, y, z, t) be the coordinates of a flat 4D space-time of signature (-I-, -|-, — , — ). 
The invariant metric reads 

ds^ = 2dxdz - 2dydt (10.1) 

The SDYM equations in 2-1-2 dimensions {e^yzt = 1) read as 

Fab = leatcdF"' (10.2) 

and are equivalently represented by 3 equations of the form 

Fte = Fy, = Fty + F,. = (10.3) 

After a dimensional reduction which is equivalent to setting 

dy = d,-d^ = 0, (10.4) 

(10.3) takes the form 

[dt -H,d^-Q] = [P, B] = 0; [H, B] = [d^ - Q, 9, - P] (10.5) 

where 

At = H; A^ = Q; Ay = -B; A, = P 

It is clear that the first equation in (10.5) is a zero curvature equation. We now choose 
the non-compact group and an embedding pattern in the form 

B=(° °); (10,6) 

« = (-.. -a) (!"■') 

where A is a constant and u = u{t,x + z). We can expand the Lie Algebra- valued fields 
H and P as 

H = H^T+ + H+T^ + ifgT-s; P = P-r+ + P+T^ + P3T-3 (10.8) 
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where t± = (ri ± ir2)/2, and ti^2,3 are the Pauh spin matrices. It is clear that the second 
equation of (10.5) gives 

P- = P3 = 0, (10.9) 

while the third equation of (10.5) gives 

H^ = -P+- Hs = -^d^{u + P'_^)-XP+ (10.10) 

where primes denote derivatives w.r.t. x. 

Finally, the first equation of (10.5) yields 3 equations 

H+ = uP+-\d^P+-]^dA{u^P+)] 9,(m + 2P+)=0; (10.11) 

ii = ^d,d,d,{u + P+) + {u-P+)d,u + 2\^P+ 

It follows that 

1 1 A 1 

P+ = -^u; H+ = --u'^ + -u^ - -u^^; (10.12) 

. _ 1 3 _ 2 

Changing the notation as 

u ^ u + -\ \ t ^ At; X + y -^ X, 
one obtains the KdV equation 

This example may be relevant towards an ultimate unification of 2D integrable models 
and 2D conformal field theories, as well as within the 4D SDYM theories which are also 
closely related to A^ + 2 strings. 

10.1 Self-Duality and Super symmetry 

Extended Supersymmetry is compatible with self-duality in 2+2 dimensions. Therefore 
the Supersymmetric self-dual Yang-Mills theory (SSDYM) is capable of generating Super- 
symmetric 2D integrable models. However a Supersymmetric generalization of the SDYM 
theory is not unique. One could either replace a gauge group by its graded version, or a 
2-1-2 dimensional space-time by superspace. 

Supersymmetric generalizations of the KdV equation in 1+1 dimensions were ob- 
tained independently by Manin and Radul [43], Mathieu [44], Bilal and Gervais [45]. 
These equations have two dynamical variables, one bosonic u{x,t), and one fermionic 
i/j{x,t), and read 

1 3 

dtu = -u^^xx + ^udxu + -{^xx)iJ (10.13) 

1 3 

<9t^ = -^ipxxx + -^dxiwp) 
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They are invariant under the A^ = 1 Supersymmetry transformations 

5u = ed.^^; 5^ = eu (10.14) 

e being a constant Grassmann parameter. Eqs.(10.13) are integrable, and can be obtained 
from the zero curvature condition associated with the graded Lie Algebra osp{2, 1) 

dtA,-d,At+[At,A4=0 (10.15) 

when the following ansatz is used for 2D Yang-Mills potentials [45]: 

2At{x,t) = \ -2u -Ux tiJx (10.16) 

\ ii/j^ iipxx + '^iuil^ / 

(0 u —iip \ 
-10 (10.17) 

i^ y 

The 2D Super KdV can be embedded into the self-duality equations by choosing the 

osp(2/l)-valued matrices H,Q as H = At{x, t), Q = Ax{x, t), and 5, P as 3 x 3 matrices 

as 



10.1^ 




P= (10.19) 

vo f^ y 

It can also be shown that the A^ = 1 and N = 2 Super KdV equations, as well as the 
A^ = 1 Super-Liouville and Super- Toda equations, can all be obtained from the N = 2 
SSDYM theory by dimensional reductions and truncations [46]. A detailed analysis is 
however outside the scope of this Article. 

11 Conclusions 

It has been our aim to present a bird's eye view of the important developments in Inte- 
grable Systems over the past few decades. What has been achieved is possibly a more 
subjective viewpoint, related to building connections between sundry topics of immediate 
interest. It has certainly not been possible to delve more deeply into the fascinating devel- 
opments in affine Toda Field Theory which seems to be a thrust area of research today. We 
refer to the excellent lecture series by Corrigan [48] on this subject. Neither is it possible 
to present an account of the interesting link between the KdV theory and Matrix models, 
encompassing thus 2D gravity; (see ref. [25] for a readable account). Supersymmetric Toda 
Field Theories have also been given the go by. They were first studied by Evans and 
HoUowood [49], as well as by Leites et al [50]. It seems to be possible to construct Toda 
Field Theories based on Lie Superalgebras with one proviso, namely, that it is necessary 
that the Lie Superalgebra admits a purely fermionic root system. This is only possible 
for the following algebras: 

A{n, n — 1); B{n, n); B{n — 1, n); D(n.n — 1); -D(n, n); andD{2, 1, a). 
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In the generic case, A^ = 1 Supersymmetric theories are obtained, which can be formulated 
in A^ = 1 superspace. There is one special case, namely, the sl{n,n — 1) theories have in 
fact N = 2 Supersymmetry; see [49]. Recently, Brink and Vasiliev [51] have proposed a 
model generalizing A]\f Toda Field Theories based on a continuous parameter, such that 
when this parameter takes on certain discrete values, the model reduces to the ordinary 
Aff Toda Theories. More recently, Wyllard [52] has worked out a WZNW reduction of 
these generalized theories, and has also attempted a Supersymmetric generalization [53] 
of the same. 

One could also picture the affine Toda theories as integrable deformations of the 
conformal Toda theories. As an example, by adding an extra simple root to the Ai 
Toda theory, one obtains the affine Toda theory, which is also the sinh-Gordon theory. 
General integrable deformations have been investigated by Zamolodchikov, among others, 
as an interesting field. Toda theories also appear in many other diverse areas of theoretical 
physics, e.g., ID discrete versions appear in the physics of monopoles [54]. Further, certain 
3D continuous Toda systems are relevant to the classification of hyper-Kahler metrics in 
4D [55]. Finally, it also appears that Toda Field Theories are relevant to M-Theory, the 
conjectured all-in-all Theory encompassing all String Theories [56]. 

I am deeply grateful to Prof. S.K. Malik, for providing me an opportunity to write 
this topical Article on the subject of Toda Field Theories. The literature on this subject 
is quite vast, but apart from my own interest in Lie Algebras and super- algebras (which 
often draws me into this area because of their obvious relevance to this subject), I have 
benefitted greatly from some concentrated literature [22,25,56,57]. I am especially grateful 
to Prof. Ashok Das for his encouragement, and for sparing time to go through the 
manuscript. I would also like to thank Jens Fjelstad for help with LaTeX. 
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